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We present an analysis on the exclusive rare radiative decay modes S;, — t- S7 
and Hfe — )• H7. The transition form factors which parameterize these processes are 



CN , calculated using QCD light-cone sum rules. The decay widths we predict are r(St — t- 

c3 : S7) = (7.21 ± 0.04) X 10-is GeV and r{Eb H7) = (1.34 ± 0.07) x 10"^^ GeV. The 

Branching ratio of — )• H7 is predicted to be Br(Hfc — )• H7) = (3.03 it 0.10) x 10~^. 

PACS numbers: 14.20.-c, 11.25.Hf, 11.55.Hx, 13.40.-f 
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I. INTRODUCTION 

Heavy flavor physics plays an important role both in the precise test of the standard 
model in the relatively high energy region and in the investigation of the hadronization 
of quarks at the low energy. Hence, a lot of effort has been paid into it and a large 
^ ! amount of experimental data have been accumulated Theoretically, much progress 

O ■ has been made in the heavy flavor meson sector for its comparatively simple structure 
2 ■ while knowledge about baryons appears to be limited. Although much literature have 
been provided to decipher these heavy flavor states (such as Refs. jsl-lsl), a deep under- 
standing of them undoubtedly demands the information on the dynamical details which 
are encoded in various decay modes jol, [l^. Among these modes, the rare radiative de- 
cay processes of the b-baryons are important in that they are not only the ways to study 
the Cabibbo-Kobayashi-Maskawa matrix elements Vts and Vtb which are closely attached 
with the dynamics inside the baryons, but also the ways to probe new physics beyond the 
standard model. 

These types of processes (such as 6 — )■ 57), which are forbidden at the tree level in 
the standard model of electroweak theory, are induced by the flavor- changing neutral cur- 
rent (FCNC) of the 6-quark. Their amplitudes are dominated by the one-loop diagrams 
with a virtual top quark and a W boson, and thus are strongly suppressed by Glashow- 
Iliopoulos-Maiani mechanism. The relative b-meson rare radiative decay modes have been 



investigated experimentally since the early 1990's lll-ll4|. while not so many experimental 



data are available for the corresponding b-baryon processes. Theoretical studies on the 



exclusive processes are available for both b- mesons and b-baryons |15l-ll9| , despite the fact 



that the dynamics of the b-baryons decays are far less clarified in comparison with that of 
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the b-mesons. However, most of the existing hterature is about the process Af, — )■ A7, the 
branching ratio of which has been predicted to be Br{Ah — > A7) < 1.3 x 10~^ experimen- 
tally [20|. Unfortunately, this decay mode is not expected to be measured easily in the 
experiments due to the fact that the final state A baryon is of neutral charged, as argued in 
Ref . j2l| . For this reason, we turn to study the possible decay modes of other octet heavy 
baryons Sf, and Sj,, in which charged final states arise and may be easily tested in experi- 
ments. It has been estimated early in the 1990's that an amount number of b-baryons may 
be produced at the cm. energy level of the LHC |22|. Thus we can expect that these rare 
decay modes could be measured by the LHC experiments in the near future, the updated 
energy of which is expected to be ~ 14 TeV. 

The remainder of this paper is organized as follows. We give an introduction to the 
exclusive rare decay mode 6 — )■ 37 and derive the formula of the decay widths in Sec. [Tll 
Then the light-cone QCD sum rules for the relative transition form factors are derived in 
Sec. mil Finally, Sec. HV] is devoted to the numerical analysis and a summary is given at 
the end of this section. 



II. PARAMETRIZATION OF THE TRANSITION FORM FACTORS 

In the standard model, the process of the exclusive rare decay b s'j can be described 
by the following effective Hamiltonian 23| : 



H,ff{b ^ .7) = -4^V*V,,CMOr{^^), (1) 

with 

Or = a^, {m,R + m,L) h (2) 

where L/R = (l=p 75)72 and F^^'^ is the field strength tensor of the photon. G^t^ is the 
Fermi coupling constant and Crlfi) is the Wilson coefficient at the scale /i. Considering 
the general form beyond the standard model, O7 can be represented as 

Or = i^rrihS a^, {gy + 75^?^) h (3) 

32 TT^ 

The decay amplitude is given by the expectation value of the effective Hamiltonian 
between the initial and final states at the hadron level 

M{X, ^ X7) = {X^\H,ff\X,) , (4) 

where X stands for the baryon involved in the process. 

By considering the Lorentz structure, the contribution of the hadronic part to the pro- 
cess, which is written as the hadronic matrix elements, is generally parameterized in terms 
of the following form factors: 

(X,(P') |j.|X(P)) = X,{P') [/i7, - hia.^q^' + hq, - {gij, + g2^(r,,q^ + gsq,h,]XiP), (5) 
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where and X are the spinors of the baryons and the weak current jj, is defined as 

ju{x) = ib{x)(7^^{l - 75)g''s(x). (6) 

In fact, form factors /s and do not contribute to the process due to the conservation of 
the vector current. Therefore, the form factors we need to calculate are fi{gi) and f2{g2), 
which can be determined from the QCD light-cone sum rules. It is noted that the processes 
are only related to the form factors at the point = 0, thus we just consider this case in 
the following analysis. 

With the form factors defined above, the decay width is represented as 

,(^^ ^ ^ Gl\V..V^^a.„^CM^^ ^M^^ ' ^^^^ ^ 

III. LIGHT-CONE SUM RULES FOR THE FORM FACTORS 

Now we apply the light-cone QCD sum rule approach to calculate the transition form 
factors fi{gi) and f^id^)- The interpolating currents to the heavy baryons are chosen as 
= e'^'1gH0)C^g^(0)]75;^fe'(0) for and js,(0) = e^^'=[s^(0)C/6^(0)]75/g'(0) for E^, 
respectively. Herein q stands for u or d quark, C is the charge conjugaion matrix, and z is 
the vector defined on the light-cone z"^ = 0. The normalization of these currents is defined 
by the parameters /s^ and f^^: 

(0|jsjs,(p')) = hs^-n^^bin, 

(Ob-^jS,(P')) = h,{z-PV^,{P'). (8) 

In the following part, we will take — )■ as an example. Our starting point for 
calculating the form factors is the correlation function 

T, = tJ rfW^(0|js,(o)J.(^)|S(P, s)) (9) 

at = and with Euclidean ml — P'^ of about several GeV^. Following the standard 
procedure of the light-cone sum rule method, we need to express the correlation function 
both phenomenologically and theoretically. By inserting a complete set of intermediate 
states and using the definitions ([5]) and ([8]), the phenomenological side is represented as 

^"^TAP, q) = ^^j^'}'^pl\ fJ - f2M - 9iH + 92Ml^]nP) + (10) 

where stands for the continuum contributions. The correlation function <^ is con- 
tracted by z'^ to remove contributions proportional to the light-cone vector z'^ which is 
subdominant on the light-cone. 
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On the other hand, the theoretical side is obtained by contracting the heavy b quarks in 
the correlation function and using the distribution amplitudes presented in Ref. 24-2^. To 
make the paper self-contained, we present in the Appendix the definition and the explicit 
expressions of the distribution amplitudes of E and S used in this paper. After assuming 
the quark-hadron duality and performing the Borel transformation, we arrive at the final 
light-cone sum rule of the form factor /2(0): 



=0 



where s = (1 — a^)M'^ + ml/as, M is the mass of the final baryon, and is the Borel 
parameter which is introduced to suppress the contributions from the higher resonances 
and the continuum states. Our calculation shows that fi = gi = and f2 = 92- In Eq- 
( ITTj) . the following abbreviations are used for convenience: 

pi— as 



B^{a3) = (2^i-^2-V3-\4-V5)(a3X 

S2(«3) = {-Vi + V2 + V3 + V^ + V^-VQ){a3). (12) 

The distribution amplitudes with tildes which come from the integration by parts in 0:3 
are defined as 

Vii^as) = / da^ / daiVi{ai,l — ai — a^ja^), 
Jo Jo 

Vi^as) = / da'^ / da'^ / da2Vi{ai,l — ai — a'^, a'^). (13) 
Jo Jo Jo 

The same procedure is also carried out to calculate the transition form factors of the 
process Eh — )■ S7. We obtain the final sum rule as follows: 

/Hj2(0)e ^ = / da2e ^ {Co{a2) + -^C^ia2) + ^-^C2{a2) - —^C3{a2)\ 
M^al^e'^ r , , 1 , , , , 

2 /l#2 I 2 1 <-^U"20j H (^2{a20) - T72"<--3ltt20j 

a^oM^C3(a2o) \ 
t2o al^M"^ + m2 

where s' = (1 — a2)M^ + m^/a2, M is the mass of H, and the following abbreviations are 
used: 

-1-02 



pi-— 012 

Co(«2) = / (iaiTi(ai, 0:2, 1 — «! — 02), 
Jo 
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CiM = {2T,-T2-T,-2Tj-2Ts){a2), 
CsM = {f2~fs-f^ + f, + f-r + fs){a2), 

CsM = i-fi + f2 + n-fe + 2fT + 2fs)ia2). (15) 
The functions with tildes are defined as 

Ti(a2) = / dao / daiTi(ai,a2,l — ai — a'2 



12 I u.^i±iyuLi,UL2, ^ "1 u;2;, 
JQ Jo 

^j(a2) = / (^(^2 \ / da2Ti{ai, 02^1 — ai — 02). (16) 

Jo Jo Jo 

IV. NUMERICAL ANALYSIS AND THE SUMMARY 

Before the numerical evaluation of the sum rules ( fTTl) and dHD, we need to determine 
the input parameters. Two important parameters are the decay constants /s^ and /s^, 
which can be calculated with the QCD sum rule approach. Using the same expressions 



in Refs. [25] and 26| with the replacements — > rrih for /^^ and rric — > mj, for /h^, we 
get the estimations /s, = (6.18 ± 0.03) x lO^^ GeV^ and /=, = (3.32 ± 0.46lx 10"=^ GeV^ 
Other input parameters needed in our calculation can be read from Ref. \2l\ : 

rrib = 4.8 GeV, = 0.15 GeV, Ms = 1.189 GeV, 
Mh = 1.314 GeV, Ms, = 5.729 GeV, M^, = 5.81 GeV, 
Vts = 0.0403, Vtb = 0.9992, (17) 

and 

aem = 1/137, Gf = 1.166364 x 10"^ GeV"^ Cjinib) = -0.31. (18) 

An important step in the numerical analysis of the QCD sum rules is to determine the 
Borel mass parameter M^ and the continuum threshold Sq. The continuum threshold Sq can 
be chosen by demanding that the continuum contribution is subdominant in comparison 
with that of the ground state which we are concerned about. Simultaneously, the resulting 
form factors should not vary drastically along with the threshold. Thus sq is generally 
connected with the first resonance which has the same quantum numbers as the particle we 
care about. Here we fix the threshold Sq in the region 39 GeV^ < < 41 GeV^. As for the 
Borel parameter M^, which is introduced to suppress the higher resonance contributions 
efficiently, we also demand that the higher twists contributions are less significant and the 
form factors should vary mildly along with it. Our calculation shows that the working 
windows can be chosen properly in the region 8 GeV^ < M| < 11 GeV^ for ^ S7 and 
9 GeV^ < M| < 12 GeV^ for S7. 



Using the distribution amplitudes given in Refs. 25| and 26|, we obtain the form factors 



at the zero momentum transfer /2(5'2)(0) as functions of the Borel parameter M^, which 
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FIG. 1: The dependence of the form factors /2(0)'s on the Borel parameter with sq 
39, 40, 41 GeV^ from the top down. 
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FIG. 2: The contributions to the form factors /2(0)'s from different twists on the Borel param- 
eter with So = 40GeV2. 



are displayed in Fig. [T] . We have also analyzed the contributions from the distributions 
of different twists, which are shown in Fig. [2J The results show that the contributions 
of the leading-order and next-to-leading-order twists are dominant while the contributions 
from higher twists are suppressed efficiently. This implies that the light-cone expansion is 
reasonable in the cases we considered in this paper. 

By using of the form factors we have estimated above, the decay widths of the processes 
can be easily evaluated with the formula which turn out to be F(Sfc — )■ S7) = (7.21 ± 
0.04) X 10-^^ GeV and F(Sfe ^ S7) = (1.34 ± 0.07}x lO^^^ GeV. Although the mean life 
of has been estimated experimentally 0, Q, 28, 29 1, here we use the average value given 
in Ref. 27|] to estimate the branching ratio of the process — t- S~7 which turns out to 
be i?r(S^ — > S~7) = (3.03 ± 0.10) x 10~^. The errors in the widths come from the choices 
of the threshold, the sum rule windows, and the uncertainties in the decay constants /sj, 
and /sj^. It is worth noting that errors from other sources are not considered here because 
the sum rule method itself brings in an amount of uncertainties (about 20%), which makes 
it less significant to take into account the errors of the input parameters. 
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TABLE I: Decay widths and Branching ratios at different points of rrih . 



mb(GeV) 


4.7 


4.8 


4.9 


r{^b S7)(x 10-18 GeV) 


6.92 ±0.03 


7.21 ±0.04 


7.26 ±0.07 


r{Eb H7)(x 10-16 GeV) 


0.98 ± 0.04 


1.34 ±0.04 


1.75 ±0.05 


Sr(H,-^H-7)(xlO-4) 


2.21 ±0.08 


3.03 ±0.10 


3.96 ±0.11 



We also investigate the sensitivity of the form factors to the variation of mb at different 
points rub = 4.7, 4.8, and 4.9 GeV. The corresponding predictions for the decay widths 
and branching ratios are given in Table [B 

In summary, we have investigated the exclusive rare decay processes S7 and 

Hfo — 7- S7. The corresponding transition form factors are estimated through the light-cone 
QCD sum rule approach and the decay widths of these processes are predicted to be TiJ^b — t- 
S7) = (7.21 ± 0.04) X 10-18 GeV and T(Eb H7) = (1.34 ± 0.04) x lO^i^GeV. We also 
estimate the branching ratio of — t- S-7, which is i?r(S^ — )■ S-7) = (3.03 ±0.10) x IQ-^. 
As we can see, our prediction is larger than the theoretical estimations for the branching 
ratio of the A^, — )■ A7. Therefore it is reasonable to assume that this mode may be tested 
easily, provided that a good source of T,b or Eb is available in the future experiments, such 
as the LHC experiments. 
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Appendix 

In the following we give the distribution amplitudes of S and S used in the paper. In 
general, the distribution amplitudes are defined by the matrix element of the three-quark 
operator as 

me''%a{aiz)q,^{a2z)q,';{asz)\X{P)) 
= V, (^C)„, (75X+)^ ± V2 (^C)„, (75X-)^ ± (7^C)^, (7^5X+)^ 

+yM (7xC)„, (7^5X-)^ ± V,— i^C)^, (75X+)^ ± —V, i^C)^^ (75X-)^ 
+Ts- {ia,,C)^p (75X+) ± T4- {ia.,C)^^ (75X-) 



8 



/Vf2 /\/f2 



72 ' \ / 7 

where M is the mass of the baryon X and C is the charge conjugation matrix. Note that 
the other Lorentz structures which do not contribute to the calculations are omitted. For 
each distribution amphtudes Fi = Vi ,Ti defined above, it can be presented as 

/—ipz Y] Xiai 
Vxe > F{xi), (20) 

with the relationship 0<a:;j<l, ^Xi = l, and Xi corresponds to the distribution of the 

i 

baryon momentum on the quarks. The integration measure is defined as 



J Vx = J dxidx2dx^8{xi + X2 + x^ — 1) . 



(21) 



The distribution amplitudes can be expanded with a conformal spin. The detailed 



process is referred to Refs. 24l426l| . The explicit expressions of the distribution amphtudes 
are collected below: 

Vi{Xi) = 12OXiX2X303, V2(Xi) = 24xiX204 , 

Vsix,) = 12x3(1 - X3)^° , V,ix,) = 3(1 - X3)V^° , 

V,{x,) = 6a;30°, V,{xi) = 2(j)l 



Ti{Xi) = 120X1X2X303°, T2(Xi) = 24xiX20!° 



4 5 



T3(X,) = 6X3(1 - X3)(e4° + C') , Ux,) = -^(Xi + X2)(e;° + ^5) 

nix,) = 6x30;°, r6(x,) = 20;° 

T,{x,) = 6x3(1 - X3)(e;° - ^4°) , Tsix,) = lix,+ X2)(e;° - e,) . m) 



The parameters in the expressions are as follows: 



0^ = <t>l = fx. r, = i^l = -{fx-Xi), 

<Pl = <Pt = lifx + Ai), 0:° = 0^° = -^5° = ^(4A3 - A2) , 

< = ^4° = ^(8A3 - 3A2), 0'° = = ^A2 , 

6 6 

= ^(12A3-5A2). (23) 

The nonperturbative parameters fx, Ai and A2 are determined with QCD sum rules to be 

/s = (9.4 ± 0.4) X 10"^ GeV^ Ai = -(2.5 ± 0.1) x 10^^ GeV^ 

A2 = (4.4 ± 0.1) X 10-2 GeV^ A3 = (2.0 ± 0.1) x 10"^ GeV^ (24) 



1 
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for E and 



h 
A2 



(9.9 ± 0.4) X 10~^ GeV^ 
(5.2 ±0.2) X 10"^ GeV^ 



Ai = -(2.8 ±0.1) X 10-2 GeV^ 
A. = (1.7 ±0.1) X 10^2 GeV^ 



(25) 



for S. 
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